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The in-plane thermal conductivity k of the iron selenide superconductor FeSex {Tc — 8.8 K) were 
measured down to 120 mK and up to 14.5 T (~ 3/4i7c2). In zero field, the residual linear term 
Ko/T at r is only about 16 fiW K^^ cm^^, less than 4% of its normal state value. Such a 
small Ko/r does not support the existence of nodes in the superconducting gap. More importantly, 
the field dependence of kq/T in FeSe^, is very similar to that in NbSe2, a typical multi-gap s-wave 
superconductor. We consider our data as strong evidence for multi-gap nodeless superconductivity in 
FeSe^;. This kind of superconducting gap structure may be generic for all Fe-based superconductors. 
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Just as CuO-plane is the basic building block of 
high-Tc cuprate superconductors, the FeAs-layer is the 
basic structure of the newly discovered FeAs-based 
high-Tc superconductorsiii^i^ii^i^i^ii The FeAs-layer con- 
sists of a Fe square planar sheet tetrahedrally coordi- 
nated by As. However, unlike the rigid CuO-plane in 
cuprates, partial substitution of Fe by Co or Ni, or 
As by P within the FeAs-layer can effectively induce 
superconductivity^i^ii^iiiii^ In this sense, the discovery 
of superconductivity in binary FeSe^; (Tc ~ 8 K) is of 
great interests, since it only contains the superconducting 
FeSe-layer which has identical structure as FcAs-laycr, 
and the Se deficiency may cause the superconductivity.^"^ 
More remarkably, the onset Tc can be enhanced to as 
high as 37 K for FeSe^; under high pressure r^^i^ii^ which 
further implies that superconductivity in FeSe^; may have 
the same mechanism as in FeAs-based superconductors. 

For this new family of high-Tc superconductors, the 
pairing symmetry of its superconducting gap is a key to 
understand the mechanism of superconductivity. Exten- 
sive experimental and theoretical work have been done 
to address this important issue for FeAs-based supercon- 
ductors (for a theoretical review, see Ref. 17; for an ex- 
perimental review, see Ref. 18). Although there is still no 
consensus, more and more evidences point to multi-gap 
nodeless superconductivity, possibly an unconventional 
paring mediated by antiferromagnetic fluctuations;^^ 
For the prototype FeSe^; superconductor, however, there 
were very few experiments to study the superconduct- 
ing gap structure. This is due to its relatively lower Tc 
and lack of sizable high-quality single crystals The 
measurements of in-plane magnetic penetration depth for 
polycrystal FeSco.ss are in favor of anisotropic s-wave su- 
perconducting gap or two gaps (s 4- s).^^ To clarify this 
important issue, more experimental work are needed for 
FeSe^; superconductor. 

Low-tempcrature thermal conductivity measurement 
is a powerful tool to study the superconducting gap 
structure;^ The field dependence of the residual thermal 



conductivity kq/T for BaNi2As2 (Tc = 0.7 K) is consis- 
tent with a dirty fully gapped superconductivity)^ For 
Bai_^K2;Fe2As2 (Tc ~ 30 K) and BaFe1.9Nio.1As2 (Tc 
— 20.3 K), a negligible kq/T was found in zero field, in- 
dicating a full superconducting gap.^^'^^ However, k{T) 
was only measured in magnetic fields up to If = 15 T 
(~ 1/4 HC2), thus can not show clearly whether the su- 
perconductivity has multi-gap character in FeAs-based 
superconductorS i^^'^^ 

In this paper, we measure the thermal conductivity k 
of a FcScj; single crystal with Tc = 8.8 K down to 120 mK 
and up to 14.5 T (~ 3/4 i?c2) to probe its superconduct- 
ing gap structure. In zero field, kq/T is about 16 
/ cm, less than 4% of its normal-state value. Such a 
small Ko/T should not come from the nodal quasiparticle 
contribution. It may simply come from the slight over- 
estimation when doing extrapolation, due to the lack of 
lower temperature data. The field-dependence of kq/T is 
very similar to that in multi-gap s-wave superconductor 
NbSe2. Based on our data, it is evident that FeSca; is a 
multi-gap nodeless superconductor. 

FeSca; single crystals with nominal formula FeSeo.82 
were grown via a vapor self-transport method. The ab- 
plane dimensions of as-grown crystals ranges from a few 
hundred /im to 1 mm. Energy Dispersive of X-ray (EDX) 
microanalysis (Hitachi S-4800) show that the actual FeiSe 
ratio is very close to 1:1 in our FeSej; single crystals. The 
nominal formula FeSeo.82 was used in the initial work by 
Hsu et al.p^ However, the actual superconducting phase 
was later determined to be FeSeo.99±o.o2 in Ref. 27 and 
FeSeo.974±o.oo5 in Ref. 28. Therefore the EDX result of 
our FeSca; single crystals is consistent with these two later 
reports. 

The ac magnetization was measured in a Quantum De- 
sign Physical Property Measurement System (PPMS). 
An as-grown single crystal with dimensions 1.0 x 0.40 
mm^ in the plane and 190 /im thickness along the c-axis 
was selected for transport study. Contacts were made 
directly on the sample surfaces with silver paint, which 



2 



o 

a 



c\i 



120 
100 
80 
60 
40 
20 

0.0 
-0.2 
-0.4 
-0.6 
-0.8 
-1.0 



0.14 








10 20 
T(K) 



30 




-^1.47(^^1.47^ 



FIG. 2; (Color online) Temperature dependence of the in- 
plane thermal conductivity for FeSex single crystal in zero 
field. The solid line represents a fit of the data to k/T = 
a + bT°-^. This gives the residual linear term kq/T = 16 ± 
2 AiW cm-\ 



FIG. 1: (Color online) (a) In-plane resistivity p(T) of FeSca; 
single crystal in H = and 14.5 T magnetic fields along the 
c-axis. The solid line is a linear fit of p{T) from 8 to 30 K, 
which gives the residual resistivity po = 57.9 p^ cm in H = 
14.5 T. (b) Normalized ac magnetization. 



were used for both resistivity and thermal conductivity 
measurements. The typical contact resistance is a few 
ohms at room temperature and 1.5 K, which is not as 
good as that on Bai_a;Ka;Fe2As2 and BaFe1.9Nio.1As2 
single crystals j^^i^^ In-plane thermal conductivity was 
measured in a dilution refrigerator using a standard one- 
heater-two-thermometer steady-state technique. Due to 
the small size of the sample and the non-ideal contacts, 
good thermalization between sample and the two RUO2 
thermometers can only be achieved down to 120 niK. 
Magnetic fields were applied along the c-axis and per- 
pendicular to the heat current. To ensure a homogeneous 
field distribution in the sample, all fields were applied at 
temperature above Tc- 

Fig. la shows the in-plane resistivity of FeSe^ single 
crystal in iJ = and 14.5 T magnetic fields. The middle 
point of the resistive transition is at Tc — 8.8 K in zero 
field. The 10-90% transition width of our crystal is as 
broad as the powder sample^i^ which has been noticed 
in Ref. 21. Above Tc, p{T) manifests a very good linear 
dependence up to 80 K, similar to the powder sampleJ^ 
A hnear fit of p{T) gives the residual resistivity po = 57.9 
pfl cm in H = 14.5 T, which is about 1/4 the value of 
powder sample.™ 

To estimate the upper critical field i?c2(0) which 
completely suppresses the resistive transition, we define 
Tc{onset) at the temperature where p{T) deviates from 



the linear dependence, and get Tc{onset) = 11.9 and 6.3 
K ioi H ~ and 14.5 T, respectively. Using the rela- 
tionship Hc2/Hc2{0) = 1 - (Tc/r,(0))2, we get ^,2(0) = 
20.1 T. Note that i7c2(0) = 16.3 T was estimated for the 
powder sample, in which Tc was defined at the middle 
point of the transition. ^'^ 

In Fig. lb, the normalized ac magnetization for FeSe^, 
single crystal is plotted. The positive ferromagnetic back- 
ground has been attributed to the existence of Fe impu- 
rity in the FeSe^; powder sample.— However, no iron, 
iron oxide, or iron silicide impurities were detected in 
our crystals;^ therefore the ferromagnetic background 
likely results from the magnetic Fe cluster promoted by 
Se vacancies!^ 

Fig. 2 shows the temperature dependence of the in- 
plane thermal conductivity for FeSea; in zero field. To 
extrapolate the residual linear term kq/T, we fit the data 
to k/T = a + 6r"^i,30.3i aT and 6T" represent 

electronic and phonon contributions, respectively. In Fig. 
2, the data from 120 mK to 0.7 K can be fitted (the solid 
line) and gives kq/T = 16 ± 2 cm~^, with a — 

2.47. 

Such a value of kq /T is slightly larger than the experi- 
mental error bar ± 5 pW cm~-'^i^ However, it is still 
fairly small, less than 4% the normal state Wiedemann- 
Franz law expectation k^q/T = Lq/pq = 0.423 mW 
cm~^, with Lq the Lorenz number 2.45 xlO^® W fl 
and Po — 57.9 pfl cm. For unconventional superconduc- 
tors with nodes in the superconducting gap, a substantial 
kq/T in zero field contributed by the nodal quasiparti- 
cles has been found.— 1^ For example, for overdoped d- 
wave cuprate superconductor T12201 with Tc = 15 K, 
Ko/T = 1.41 mW K-2 cm-i, about 36% kno/TM For 
p-wave superconductor Sr2Ru04 with Tc = 1.5 K, kq/T 
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FIG. 3: (Color online) Low-temperature thermal conductivity 
of FeSe^, in magnetic fields applied along the c-axis {H = 0, 
f, 4, 9, and 14.5 T). The solid lines are n/T = a + bT"'^ fits. 
For // = 4, 9, and 14.5 T, the electronic contribution becomes 
more and more dominant and the data get less smooth, there- 
fore a is fixed to 2.47. The dashed line is the normal state 
Wiedemann- Franz law expectation at T — > 0, namely Lo/po, 
with Lo the Lorenz number 2.45 xlO"^ W Q K"^. 



= 17 mW K^^ cm^-'^, more than 9% knq/T for the best 
samplci'^'^ We also note that kq/T in zero field are all 
negligible in closely related superconductors BaNi2As2, 
Bai_:rK:rFe2As2, and BaFei.9Nio.iAs2.^'25^ Therefore, 
it is unlikely that kq/T = 16 ± 2 R-^ cm"! in FeSe^^ 
single crystal conies from the nodal quasiparticles. Since 
no impurity phases were detected, such a small hq/T may 
simply come from the slight overestimation when doing 
extrapolation, due to the lack of experimental data below 
120 mK. 

Below we turn to the field dependence of kq/T in 
FeScj;. Fig. 3 shows the low-temperature thermal con- 
ductivity of FeScj, in magnetic fields applied along the 
c-axis {H = 0, 1, 4, 9, and 14.5 T). For = 1 T, the 
data is also fitted to k/T = a + bT"'^, and gives kq/T 
= 47 ± 2 K-2 cm-\ with a = 2.47. For H = 4, 
9, and 14.5 T, the electronic contribution becomes more 
and more dominant and the data get less smooth, there- 
fore a is fixed to 2.47 in the phonon term feT". From Fig. 
3, even higher magnetic field is needed to increase k/T 
to its normal state value. 

In Fig. 4, we put the normalized k,q/T{H) of FeSe^ 
together with the clean s-wave superconductor Nb,-^"* the 
dirty s-wave superconducting alloy InBi,'^^ the multi- 
band s-wave superconductor NbSe2'^® and an overdoped 
sample of the d-wave superconductor Tl-2201,'^^ plotted 
as a function of H/Hc2- For a clean (like Nb) or dirty (like 
InBi) type-II s-wave superconductor with isotropic gap, 
Kq/T should grow exponentially with field (above -ffci)- 
This usually gives negligible kq/T for field lower than 
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FIG. 4: (Color online) Normalized residual linear term kq/T 
of FeSea; plotted as a function of H/ Hc2- Similar data of the 
clean s-wave superconductor Nb,'''* the dirty s-wave supercon- 
ducting alloy InBi,-^ the multi-band s-wave superconductor 
NbSe2,— and an overdoped sample of the d-wave supercon- 
ductor T1-2201S are also shown for comparison. 



i?c2/4. For the d-wave superconductor Tl-2201, kq/T 
increases roughly proportional to ^/H at low field due 
to the Volovik effect. ■^^ By contrast, for multi-gap super- 



conductors NbSe2 and MgB; 



36.38 



magnetic field will first 



suppress the superconductivity on the Fermi surface with 
smaller gap, and give distinct shape of kq/T{H) curve, 
as seen in Fig. 4. 

From Fig. 4, the kq/T{H) of FeSca, manifests almost 
identical behavior as that of multi-gap s-wave supercon- 
ductor NbSe2. For NbSe2, the shape of kq/T{H) has 
been quantitatively explained by multiband supercon- 
ductivity, whereby the gap on the F band is approx- 
imately one third of the gap on the other two Fermi 
surfacesi^ Therefore, we consider our data as strong evi- 
dence for multi-gap nodeless superconductivity in FeSe^: . 
Note that in the two-gap s -I- s-wave model to describe 
the in-plane penetration depth data, the magnitude of 
the two gaps are 1.60 and 0.38 meV, respectively^-^?' The 
ratio between these two gaps is about 4, close to that 
in NbSe2, thus supports the multi-gap scenario from our 
thermal conductivity results. 

So far, there is still no experiment to directly mea- 
sure the superconducting gap in Fei+j^Tei_2:Se2: sys- 
tem. Density functional calculations show that the 
electronic band structure of FeS, FeSe, and FeTe are 
very similar to the FeAs-based superconductors.— In 
doped BaFe2As2, multi-gap nodeless superconductivity 
has been clearly demonstrated by angle-resolved photoe- 
mission spectroscopy (ARPES) experiments . For 
hole-doped Bao.6Ko.4Fe2As2 {Tc = 37 K), the average gap 
values A(0) for the two hole pockets {a and (3) are 12.5 
and 5.5 meV, respectively, while for the electron (7 and 
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S) pockets, the gap value is similar, about 12.5 meVj ^^i*^^ 
For electron-doped BaFei.85Coo.i5As2 (Tj, = 25.5 K), 
the average gap values A(0) of hole {(3) and electron 
(7 and 5) pockets are 6.6 and 5.0 meV, respectively^^. 
The ratio between the large and small gaps is 2.3 for 
Bao.6Ko.4Fe2As2. This may explain the rapid increase 
of Ko/T{H) at low field in Bai_a;K^Fe2As2,^^ although 
magnetic field was only applied up to 1 /4 Hc2 thus could 
not see clear multi-gap character as in our FeSca; single 
crystal. 

In summary, we have measured the low-temperature 
thermal conductivity of iron selenide superconductor 
FeScj: to investigate its superconducting gap structure. 
A fairly small kq/T at zero field and the dramatic field 
dependence of kq/T give strong evidence for multi-gap 
nodeless superconductivity in FeSe^;. Such a gap struc- 



ture may be generic for all Fe-based superconductors. 
More experiments are needed to distinguish conventional 
s-wave from the unconventional s^-wave superconduc- 
tivity in this new family of high- Tc superconductors. 
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